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Abstract 

The first order post Newtonian scheme in multiple systems presented by Damour-Soffel-Xu is extended 
to the second order one for light propagation without changing the advantage of the scheme on the linear 
partial differential equations of potential and vector potential. The spatial components of the metric tensor 
are extended to the second order level both in the global coordinates (c/ij/c 4 term) and in a local coordinates 
(Qab/c 4 term). The equations of qij (or Q a b) are deduced from Einstein field equations. The linear relation- 
ship between qtj and Q a t are presented also. The 2PN equations of light ray based on the extended scheme 
are deduced by means of the iterative method. We also use parametrized second post Newtonian metric 
tensor to substitute into the null geodetic equations to obtain the parametrized second order equations of 
light ray which might be useful in the observation and measurement in the future space missions. 
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1 Introduction 



Recently, in terms of advanced technology a series of space missions are proposed and planned to launch, e.g. 
LISA (Laser Interferometer Space Antenna) pQ, aimed to detect low-frequency (10 -4 -l Hz) gravitational wave 
with a strain sensitivity of 4 x lCP 21 /(-ffz) 1//2 , GAIA (Global Astrometric Interferometer for Astrophysics) [2] 
planned to measure the position and the parallaxes of celestial bodies with the precision of few /xas (micro arc 
second) , ASTROD ^ and ASTROD-I 0] (Astro-dynamical Space Test of Relativity using Optical Devices) in 
a lower frequency range and higher sensitivity for the gravitational wave detection compared to LISA, LATOR 
(Laser Astrometric Test Of Relativity) to measure the bending of light near the sun to an accuracy of 0.02 /xas 
5 and so on. In such a requirement of experimental precision the second order post Newtonian contribution 
to light propagation has to be taken into account. On other hand, when the light pass nearby the sun or any 
planet, the potential produced by 2 N multipole moment is V mono (l/r) N [H], where V mono is the potential created 
by the mass monopole, / the length deviated to monopole and r the distance of the action. For the quadrupolc 
of the sun (N = 2) (l/r) 2 < e 2 , where e 2 is the small parameter of post Newtonian expansion (nearby the 
sun e 2 ~ 10~ 6 ), I can be estimated from the oblateness of the sun [7J- Therefore if we consider the second 
order post Newtonian (2PN) problem, the hrst post Newtonian (1PN) quadrupole has to be taken into account. 
Because the trajectory of light ray normally is calculated in the global coordinates, but the relativistic multipole 
moments have to be computed in the local coordinates. In 1PN level, the problem has been thoroughly solved in 
the scheme presented by Damour-Soffel-Xu |8l l9l 1101 111] , in which a complete 1PN general relativistic celestial 
mechanics for N arbitrarily composed and shaped, rotating deformable bodies is described. Their scheme is 
widely abbreviated as the DSX scheme [T2]. A similar results has been deduced by means of a different way: a 
surface integral derivation |13| . 

As for light propagation, we should also mention the other methods. The method of Lorentz covariant 
theory employed to study the problems related to the propagation of light ray in a gravitational field is based 
on the solution of the null geodesic equations by means of the first Minkowskian approximation ^] E ■ In 
the second order approximation the Lorentz covariant theory has to be extended also. Numerical simulation is 
also an important tool to study the light propagation in the gravitational field of the moving bodies |16| . A 
more rigorous approaches (e.g. the second post Newtonian (2PN) scheme) might be helpful for such kind of 
simulations. 

The 2PN contribution for light ray has been discussed for a long time (early 80 's) by Epstein and Shapiro 
|17|. Richter and Matzner |18|. and by others (later). But all of them consider only in one global coordinate 
system, therefore they can not calculate the relativistic contribution from multipole moments which should be 
calculated in the local coordinates. As we know the relativistic theory of reference system established only after 
1991 1 10 ', although it is only 1PN approximation. 

In this paper we will extend the 1PN DSX scheme to 2PN for the discussion on the propagation of light 
ray, i.e. we will extend the metric in global coordinates (and G a b in local coordinates) to 0(6) and deduce 
the corresponding equations from Einstein field equation. In the section 3, we deduce the second order post 
Newtonian equation of light ray by means of the iterative method. In the section 4, we discuss the parametrized 
second post Newtonian formalism (PP 2 N) which is a tentatively expansion of Ref-HHj and the corresponding 
equations of 1PN light ray. Some conclusion remark is made in Sec. 5. 

2 The extension of DSX scheme 

In this section, we will extend 1PN DSX scheme to 2PN scheme for discussion on the second order contribution 
in light propagation. Our symbols and signature follow the DSX scheme. Here we summarize the notation in 
DSX paper which we will use in this paper. The signature -H — | — h is taken; spacetime indices go from to 3 and 
denoted by Greek indices, while space indices (1 to 3) are denoted by Latin indices. We use Einstein's summation 
convection for both types of indices, whatever the position of repeated indices. The flat metric is denoted by 
/ M „, with components diag (-1,+1.+1,+1) in Lorentzian coordinates. In post Newtonian expansions we shall 
often abbreviate the order symbol 0(c~ n ) simply by 0(n). The "global" ( or "common view") coordinates 
used for describing the overall dynamics of the system will be denoted by (x^) = (ct,x l ). By contrast, each of 
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the "local" coordinate systems A, used for describing the internal dynamics of each body, will be denoted by 
(X%) = (cTa,x a ). We distinguish the second part of the Latin alphabet (i,j,k, • ■ •) for global space coordinates 
from the first part of the Latin alphabet (a,b,c, • • •) for local space coordinates as done in the DSX scheme. 
A spatial multi-index containing I indices is simply denoted L (and K for k indices, etc.), i.e., L = ixi% ■ ■ - h- 



A multi-summation is always understood for repeated multi-index SlTl — Ei 1 ...j 1 Sj. 



T- 



vector, v l , its lih tensorial power is denoted by v 



Also, d L = d h . . . di, 



Given a spatial 
The symmetric and 



trace-free (STF) part of a spatial tensor will be denoted by angular brackets (or by a caret when no ambiguity 



arises): STF;, 



,A T n. 



T, 



(*!...<! ) 



-Tl- The coordinate transformation between the local coordinate system 



Xj( and the global one x^ reads (omitting a labelling index A on all quantities pertaining to the local frame): 



where 



= F(X a ) = z^(X°) + eZ(X°)Y a (X°, X b ) + C 



Y a (X°,X b ) =X a + (l/c 2 )[(l/2)A a (A 2 ) - X a (A b X b )]. 



Moreover, one has the definitions 



A a (S) = f MV e£(S) 



d 2 z» 
dr 2 



dz l 
Z dz° 



dz l 



eUS) = el— + 0(4) 



dS 



and 



e°(S) = 1 + 3 



1 /l 



r 



i 



TTT„.2 



e° a (S) 



Rli- 



jV 










4 -l 




H 










I" 













•0(5), 



0(6), 



1 

c z 



S ij + r-svV 
2cr 



K + 0(4) , 



(1) 

(2) 
(3) 
(4) 



where R l a (S) is a slowly changing rotation matrix R a R 3 a — , R a R l b = & at \ cAR l a /dS = 0(2). All of above three 
equations are already shown in Eq.(5.21b)-(5.21d) of Ref.jHj, as for Eq.Q will be extended to 0(6) later (see 
Eq.||£IJ). Now, in this paper, a fixed-star coordinates has been chosen, therefore R l a = S a , z^ represents each 
central world line £a(A° = 0). w is the external potential of body A, v l is the coordinate three- velocity of the 
central world line La measured in the global coordinate system, and V a defined by v l — R a V a or V a = Rfv 1 . 

As we mentioned early, DSX scheme is the first post Newtonian approximation for particle motion. The 
metric tensor in a global coordinates is written in the form 



goo = -cxp T + 0(6) , 



2w 
7i 



g i3 = Sij exp 
and they satisfy the conformal isotropic condition 

ffOOffij = Si. 



0(5), 



0(4). 



(5) 
(6) 
(7) 

(8) 



The metric tensor in the local coordinates A has a similar form as in global, one only needs to change the small 
letters into capital letters 



Goo 

Goa 
G a b 



— exp 

4W a 
c 3 

Sab exp 



2W 



0(5), 



/ 2W 



0(6), 



0(4), 



(9) 
(10) 
(11) 
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and they also satisfy the conformal isotropic condition 

G 0Q G ab = -S ab + 0(4) 



(12) 



The potential W and vector potential W a for body A can be divided into two parts (self part W + , and 
external part W, W a ) 



W = W+ + W, 
W a = W++W a . 



(13) 
(14) 



As in DSX scheme, we use W a to represent W (a — 0) and W a (a = a). From Einstein field equation, we 
obtain equations to be satisfied by W and W a (compare Eq.(4.3a), (4.3b) and (4.4) of Ref. 8 ) 



D# + — d T (d T W + d b W b ) = -47rGE + 0(4) , 
& 

A x W a - d a (d T W + d h Wb) = -47rG£ a + 0(2) , 
where Ox = Ax — cr 2 &^ with = d 2 /dX a dX a , and where the source terms 

t^OO _i_ rpaa rpQa 

S a = (E,E a )E 



(15) 
(16) 



(17) 



are now defined by components of the stress-energy tensor in the X a coordinate system. 
In harmonic gauge 

d T W + d a W a = . (18) 

We have (see Eq.(4.51) of Ref.®) 

□xW+ = -4^G£ Q . (19) 
By means of Blanchet-Damour multipole moments |20| we obtain the solutions 

W +A (T,X) = Gj2^-9 L [R- 1 M A (T±R/c)] + \d T A A + 0{4), (20) 



I 

l>0 

l 



l\ 

i>i 



<>l i ( /r^M^rCTiE/c)] +^e o6c a feL _ 1 [ J R- 1 5 c l_ 1 (T±fl/c)] 



-±d a A A + Q(2), (21) 



where 



aa = ^ ji+l^Xl 9 ^ 1 ^ ± r/c)] > (22) 

j3 v v bL-s^b 



l*t{T)= / «I 4i E s (T,X). (23) 



A 



In Eas. l|2U|l - H23|l the ± sign in a function of the local time, T, divided by the coordinate distance to the origin 
of the local system, R = |X|, denotes the average 



F(T ± R/c) _ 1 f F(T-R/c) F(T + R/c) 
R = 2 I R + R 



(24) 



which is the well-known time-symmetric solution, with spherical symmetry, of the wave equation. 

The "mass" (M A ) and "spin" (S A ) multipole moments of body A appearing in Ecis.l|2T}jl. J2IJ are the 
STF Cartesian tensors defined by the same expressions of the matter distribution variables [2]] for closed 
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gravitationally self-interacting system, but now restricted to an integration over the volume of body A, using 
the local-system matter variables Ea.i|17[l: i.e., 



4(2/ + 1) 



(1 + l)(2Z + 3)c 2 dT 

A 



J ftrVjT.XjJ +O(4)(Z>0), (25) 



S£(T) = / d 6 Xe ab<Cl X ^ L>a Z"{T,X) + 0(2){l > 1) . (26) 

J A 

In other hand, from the definition of the external part in the global coordinates, we have in the global coordinates 
for N bodies system 

N 

<=£ < , (27) 

where u>^f is the four potentials of the body B . From now on we will omit the labelling index A on all quantities 
pertaining to the local frame. Based on w^, we can calculate W a with the following equations 

/ 1V 2 \ 4 r 2 

w =(l + —jw + -V a W a + - \n(A° A° - A° a A° a ) + 0(4) , (28) 

vh = v'W + KW a + jiA^Ai - A\A\) + 0(2) . (29) 

Then from some calculation we can get W and W a also. Up to now we have briefly reviewed the main story 
of DSX scheme. But as we mentioned early that DSX scheme is 1PN approximation method (for the equations 
of motion of bodies). If we want to extend DSX scheme to apply to the second order post Newtonian (2PN) 
approximation of light propagation, we have to extend the metric <?y to 0(6) level rather than 0(4) |21| . i.e. 

Goo = -expf-— J +0(6), (30) 

G 0a = -^ + 0(5), (31) 

G ab = 5 ab exp (jpj +^ + 0(6) . (32) 



The contravariant metric tensor reads 



/ 2W 



G uu = -exp ^— J +0(6), (33) 

G 0a = -1^+0(5), (34) 

G ab = Sab exp(-?p) -^ + 0(6). (35) 

The metric tensor in the global coordinate system has a similar form, just change G a ^, W, W a and Q ab by g^ v , 
W, Wi and qij. 

The relations between w , Wi , and W , W a , Q ab will be studied in the following. In DSX scheme, since 
the metric in the global coordinate system as well as in every local coordinate systems has the similar form, so 
that, if we deduce any equation in one local coordinate system, then we have it in every local coordinates as 
well as in the global coordinates. Now in this paper, the situation is similar, Therefore we need only to treat in 
a local coordinate system A, but we omit the index A always. The spatial conformal isotropic condition Ea. (|12jl 
is revised as 

G ab G o = S ab - % + 0(6) . (36) 



5 



If we attribute Q a b/c 4 to 0(4), it is just the form in DSX scheme. Therefore Ea. l|32|l is an extension of Eq.JJJ. 
Qab also can be taken as a spatial anisotropic contribution in second order (see Eq.(J2EJl). 

The metric tensor Go a is taken as 0(3) here, in fact it is only as small as 0(4) because of the relatively slow 
rotation rate in the sun |181 \2'2\ . To keep the uniformity with DSX scheme therefore we take such formula of 
Goa as well as for g^. 

We have to calculate the equations satisfied by W, W a and Q a b from the Einstein field equations. On the 
first step, the Christoffel symbols can be calculated from the metric tensor 

r^o = -^ + o(5), 

C = -^ + 0(6). 

c c 1 c 

r° ab = sJ^ + ^w M + o(S), 

K„ = W [aM + -^-8 ab + 0(5) , 

T a bc = ^(5 ab W, c + 5 ac W, b -6 bc W, a ) + ^(Q ab , c + Q ac , b -Q bc , a )+0(6), (37) 

where two indices enclosed in a parentheses is denoted as symmetrization, and in a square bracket means 
antisymmetrization. 

Then we can deduce the Ricci tensor 

i? 00 = -^^-^(3dtd t W + 4dtd d W d ) + 0(6), (38) 

R 0a = {V 2 W a - d a d d W d - dtd a W) +0(5), (39) 

R ab = -^6 ab V 2 W + \ [4S ab WV 2 W + 5 ab W, u + 4W (0i 6 )it - 2W, a W, b 



+ -^{Qad,bd + Qbd.ad ~ Qab,dd ~ Qdd.ab) 



+ 0(6) , (40) 



and the scalar curvature 



R = -lv 2 V7 + ±-(WV 2 W + 6W, tt + 8W a , at - 2W, a W, a + Q ab . ab - Q aaM ) + 0(6) . (41) 
c c 1 

From Einstein field equation we can derive equations of W, W a and Q a b 

V 2 W + \(3W M + Ad t d a W a ) = -4^G£ + 0(4) , (42) 
cr 

V 2 W a - d a dbWb - d t d a W = -4^G£ a + 0(2) , (43) 
-25 ab W }U + 4(W (0l6)lt - S ab W {d ,d).t) ~ 2W, a W, b + 8 ab W A W A 

+ -^{Qad,bd + Qbd.ad — QabAd — Qdd,ab — S a bQdc,dc + S a bQdd,cc) = 8nGT ab , (44) 

where £ = (T 00 + T aa )/c 2 and S a = T 0a /c (the same definition in DSX scheme). Eq.lO and @!| are linear 
PDE, which can be solved in certain suitable gauge conditions, then Ea. H44|) will be solved also. We shall point 
out that, T ab can be expressed by £ and E a (to see Eq.(A13) of Ref.|23|). Therefore it is self-consistent within 
the framework of DSX scheme. Similar field equations in the global coordinate system can be obtained easily 
by substituting capital letters by small letters. 

From Eas. (|42() and (|43|l we see that Q a b does not appear in the field equations of W and W a , and they are 
the same as ones in DSX scheme (keeping linear equations) . The solutions of W and W a related to relativistic 
multipole moments are still valid as before in DSX scheme. 



G 



Now, we shall discuss the transformation relations of the potential, the vector potential and Q a t (or q^) 
from the global coordinate system to a local coordinate system and vice versa. The coordinate transformation 
law reads 

9 dx«dxn ■ [ ' 



Considering g , we obtain 



/ 1V 2 \ 4 r 2 

w =(l + —) W + ^V a W a + - ln(A°A° - A^ a ) + 0(4) , (46) 

where A 1 ^ = dx^/dX a , V a = Rfv 1 (i?" is a slowly changing rotation matrix and v % is the coordinate three- 
velocity of the central world line measured in the global coordinate system) . 
For g° l , it turns out to be 

m = V *W + KW a + ^(AlAl - A° a Al) + 0(2) . (47) 

In Eas. (|46() and l|47[) . Qij does not enter these equations, and the relations between w, Wi and W, W a are the 
same as in DSX scheme (to keep linear relation). As we know, all of the discussion about the theory of reference 
systems are based on the linear PDE of potentials fEa. H42|) . tylty ) and linear relation between potentials in the 
global coordinates and in the local coordinates f )46|l and H47fl ). Since these equations are the same as before (in 
DSX scheme), the theory of reference system in DSX scheme (for potential and vector potential) is valid in this 
paper. Also in Eas. l|4^|) . (|4*3|l . il4l)|) and lH7|l Q a b does not appear, therefore self-parts W + , and external 
parts W, W a can be obtained just like in the DSX scheme [HUH- 

Finally, for g % i , we get 

qij = -2W(2V 2 5 VJ - B} a RlV a V b ) - &V a WJ %J + 8i#i# W a V b + R a R{Q a b 

+2Wc 2 {A\Ai - Sa) + ^{AiAi - A\A{ + 8 VJ [1 - ln(A° A° - A° a A° a )}} + 0(2) . (48) 

Here we should emphasize that Ea. l|48|l has only formally given for completeness. 

We should point that in Ref.|23! our calculation is incomplete since in Ea. (|48|) A^ must be calculated up to 
0(6) level, but in fact A ] a is related with e J a as 

Ai = < + ^^ ( 49 ) 

where e 3 a is the value of A J a at the central world line of body A. But precision of e 3 a is only 0(4) level (see 
Eq. (£]}), therefore we need to extend e J a to 0(6). We have the formula (see Eq.(5.19) of Ref.(H]) 

0a/3 e a e b = 5 ab ■ (50) 

From Ea. (|5*U)l and (0J, we have 

<=K-^\A + ^v\V a + ^+0(6), (51) 

where 

K = \ w 2 \ A R a + + I vJ\ A V a v< - \ q tJ \ A Rl . (52) 

where \ A denote the value at the central world line of body A. In Ea. l44f) . W and W a can be taken as known 
functions. By means of harmonic gauge (d a W a + d t W = 0), Ea. (|4"4l becomes 

25 ab W,tt + 4W Mit - 2W, a W, b + SabW^W.d 

+ 2 (Qad,bd + Qbd,ad + SabQdd^cc ~ QabAd ~ Qdd,ab ~ &abQdc,dc) = 8^GT ab , . (53) 
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We also divide Q ab into self-part and external part: 



Qab = Q + ab + Qab ■ (54) 

In DSX scheme, they took the "weak effacement of post Newtonian external gravitational potentials in the local 
frame" , i.e. the "external" PN potentials W (Ta, Xa) in the local A frame vanish for all Ta times, at the origin 
of the frame (at the central world line) : 

V Ta, Wt(T A ,0,0,0) = (55) 
(see Eq.(5.12) of Ref.[5]). We assume "weak effacement" condition to be valid also for Q abl i.e. 

\/T A , Qf 6 (T A ,0,0,0) = 0. (56) 
With the weak effacement condition, we directly write out equations satisfied by Q^ b and Q ab from Eq.JHHJ): 
28 ah W% + 4W+ b) t - 2W+W+ - 2(W+W, b + W, a W+) + 6 ab W+W+ + 26 ab W+W, d 



2 ( Q ° d - 



bd ~ ^bd,ad 



' ^Q + dd. cc ~ Qlbdd ~ Q + dd.ab - Sa*Q±dc ) = ^ GTa " - (57) 



28 ab W, tt + W {a , b)tt - 2W, a W, b + S ab W, d W,d 

+ 2 (Qad,bd + Qbd,ad + $abQdd,cc ~ Qab,dd ~ Qdd,ab ~ ^abQdc,dc) = • (58) 

The weak effacement condition is automatically satisfied for Eq.JSSJ. In principle from Eqs.JSTJ and l|58(l and 
boundary condition we could solve Q^ b and Q ab numerically. But for the problem of the light propagation 
in the solar system it becomes much simple. Since the velocity of the relative motion is small inside the sun 
(v 2 /c 2 < e 2 ), especially the shape of the sun is close to a monopole, therefore if the origin of the coordinate 
system is taken as the center of the solar mass (dipole always equal to zero) , the quadrupole terms (and higher 
multipolc moments) in 2PN contribution on Q ab can be ignored, i.e. 

Qab - S ab Q/3 (59) 

in the local coordinates of the sun. As for other local coordinates of each planets 2PN contribution is too small 
(10 -19 — 10~ 18 ) to be calculated. In the local coordinates of the sun the external potential is negligible, then 
W = W+. 

Substituting Eq.JSHJ into Eq.Q and taking PN gauge (3d t W + Ad a W a = 0(2)), Q satisfies the equation 
as following: 

3 2 



.,/</ 



W z = 8ttG T aa - — - . (60) 



The right hand side of above equation has a compact support source and the solution of Eq. (|60[1 is 

Q = _ r6G(T«*-W md3x _3 w2 
J r 2 

From Ea. l|61|l it is clear that Q itself does not have compact support neither is linear in W (because of the 
last term). We should point that if we use PN gauge, W and W a still can use B-D moments to expand , but 
one needs to do the gauge transformation from the harmonic gauge to the PN gauge by means of choosing a 
suitable A in Eqs.lO and 

In the global coordinates of the solar system, we have a similar equation and solution. Since there is 
very small difference between the global coordinates of the solar system and the local coordinates of the sun 
(additional 10~ 3 in c~ 4 level) we can always ignore the distinction between q and Q of the sun. Therefore, for 
the situation of the solar system, it may be unnecessary to solve Eas. (|57(l and (|58J) directly, but in the following 
discussion we will still keep in the equations for the generality. 
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3 2PN Equations of Light Ray 



In this section we will deduce the 2PN equations of light ray in the global coordinates by the iterative method. 
Certainly, the result can be also used in the local coordinates by substituting the quantities in the global 
coordinates for the corresponding quantities in the local coordinates. For example, by replacing w, Wi, Qij, g^u 
by W, W a , Q a b, Gab- We start the basic equations of light ray: 



dx* dx v 
d 2 x^ „ dx v dx a 

l = o 

d\ 2 va d\ dx ■ 



(62) 
(63) 



where A is an "affine" parameter. Normally we replace A by time t in terms of [i = component of Eq. l|63[) (see 
Ref.0), then Eq.© and become 



dt 2 



l r o <^_ T i \ dx^dx^ 
c va dt va ) dt dt 



and 



dx* dx" 
~dT~dT 



0. 



(64) 



(65) 



In fact Eas. l|64|l and l|65|) are our real basic equations of light ray. For 1PN equation of light ray, we use 1PN 
metric (neglect g oi and all 0(4) terms in and <7 o) and corresponding Christoffel symbols Ea.lpT7)l: 



5oo 



-1 

5ij 



— + o(4) 



2uj 



o(4). 



Then Eas.ljB^I and j^St become 



2w - 1 + — = . 



d 2 x l 
Itt 2 



The solution of equations in 1PN level is 




4 dx 1 
^~dT 



dx „ 
— • Vw 

dt 



x = x + cn(i - t ) + xip 



(66) 
(67) 

(68) 
(69) 

(70) 



where Xi p is 1PN revision based on the Newtonian or the zeroth order solution (x = xq + cn(t — to)) and n is 
the directional unit vector which is given in the initial conditions. Finally, we have 



cn 



dxu 
dt 



d 2 xip 



= -2w. 



2 [Vu> - 2n (n • Vw)] 



(71) 
(72) 



Eqs.JJlJ and H72|l agree to the parametrized 1PN equations of light ray deduced by will 0, when we put 7=1 
and replace Will's U by w in DSX scheme. Based on 1PN solution of light ray (xip as a known function), we 
deduce the 2PN equations of light ray in the extended DSX scheme by using the iterative method. Substituting 
the 2PN metric tensor in global coordinates (Eqs . (J2DE2J replaced G a p, W, W a and Q a b by g^, w, Wi and qij) 
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into Eas.l(M| and ffiffy. we have 

- 2w + 
w, 1 



= c 2 
d 2 x l 



2w 2 



8wi dx l 
c 2 dt 



1 + — 

c z 



2w 2w 2 





dx 




~dl 



dt 2 



Aw 

72" 



1 



rfx 
~dt 



c^ c 

dx 1 I 3w t t 2wj dx 3 2w(j t fy dx 3 dx k w 



dt \ 2c 2 



dt 



dt dt 



t 

2 C 4 



qij dx % dx 3 
~c T ~dT~dt 



dx 



(73) 



dt 



Awi t t & w [i,j] dx 3 qij t k dx 3 dx k qjh,i dx 3 dx k 
c 2 + ~~c 2 dt c^~dT~dt ~ ~2?~dt~dT 



(74) 



We assume that the solution is x = xjv + x±p + X2P, where xjy is Newtonian solution (straight line), xip and 
X2P are 1PN and 2PN post Newtonian revision respectively. Then we have a series of relations: 



x l + cn 1 (t - t ) + x\„ + xL , 



dx i 
~dt ' 
d 2 x % 

HF 

dx 
It 

dx 1 dx 3 
1t~dT 



dx 



lp 



dXn 



2p 



dt 



dt 



d x\ ~, d 1 



2/j 



dt 2 



dt 2 ' 



2cn- 



r/x 



lp 



dt 



2cn 



dx2p 
dt 

q - - - dx 1 _ ■ dx 
c 2 n l n 3 + cn 1 —^- + en? — 
dt dt 



dx 



dt 



— + cn 



dt 



cn 



? dx\ p ^ dx\ p dx\ p 



dt 



dt dt ' 



(75) 
(76) 

(77) 

(78) 

(79) 



where we have neglected all of the terms higher than 2PN in Eqs.JTSJ) and (|79|l . Substituting Eqs.JTHJ), (|78|l and 
(|79")) into Eq.JZBJ, an d considering Ea. l(TT|l . we obtain 



rfx2p Aw 
dt ~ 



% i 

n 



1 

2c" 



dx\p 



dt 



Aw 2 



(80) 



where = quj\ + Sijq/3, and q^ is a trace- free symmetric tensor. In the global coordinates of the solar 
system, quj\ — 0, therefore q^ = Sijq/3. In Eq.JSOJ Xip can be solved from Eas. (|71|l and (|72|l . then Eg. (|80f) 
becomes 



dx 2P _ Aw l 
dt ~ ~ 



1 

2c 



dx. 



ip 



dt 



Aw 2 



q 

6c 3 



(81) 



Substituting Eqs.jnjl, C3, GHJ and int o Eq.lf73Jl. we have 



b 2P 



dt 2 



i i k ° 7 
= — W,t n H — w (j,k) n n n H — w [j,fc] n 2 WW - i 



dx u 

dt 



■ Vu> n 



4 . _ , dx, P 



1 qij,kn 3 n k + -^q jkti n 3 n k . 



(82) 



If 3y = Sijq/3 is taken, we obtain 



"2P 



A 2 



= — H — W(j^n l n 3 n k H — w^- ^n^ ^ww^., 



dx i I 

(n ■ Vttf) — ^ - —j (n • Vq) n 4 + — <?, 

(XE 3c z 



1 

6c 2 



dxu 



(83) 



Eas. l|8U|l and l|82|l are just the 2PN equations of light ray that we expected. By means of Eqs.JSOJ and Q82)). 
based on Eqs.(7TJ and J72J), we can discuss light propagation at the 2PN level by means of the iterative method. 
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In the special case of the solar system, we may use a much simple formulae (Eqs.JHU an d ffity). where Q (or q) 
is presented in Ea. (|6T)> |22|. 

The 1PN term x%p in the equations can be obtained from the solutions of the 1PN equations of light ray, 
i.e. from Eqs.JTlJ and (|72|l . Then the 2PN equations are solvable. Furthermore, if we knew ever higher order 
metric tensor (higher than 2PN), we also could get higher order equations of light ray by means of such iterative 
method. Besides, the iterative method can also be discussed on the parametrized 2PN equations of light ray, 
then we can discuss the 2PN equations of light ray in alternative gravitational theories which we will discuss in 
the following section. 



4 Parametrized 2PN Equations of Light Ray 

The extension may also be used in parametrized post Newtonian (PPN) formalism, which would be found a 
widely use for the relativistic experiments and observation |7 . In terms of more exact consideration |19|. PPN 
extension for DSX scheme (1PN) is written as 

.900 = -l + ^~^pw 2 + 0(6), (84) 
cr c 4 

5o, = -^-p^/ + 0(5), (85) 

Qij = Sa (l + ^v?j + 0(4) . (86) 

Comparing to a known parametrized second post Newtonian formalism |25j our extension for parametrized 
second post Newtonian formalism is taking a tentative form in the extended 2PN DSX scheme 

ffoo = -l + ^--^0w 2 +0(6), (87) 
c c 

5o, = -^±^ + 0(5), (88) 

1 + + "-f (2- 2 + |) + + 0(6) , (89) 



9ij — $ij 

When 7 = /3 = e = 7/=l, then the metric tensor return to the case in general relativity. If we ignore all of 1/c 4 
terms in , our results agree with the one of Ref . 

From the parametrized 2PN metric tensor we can get its contravariant tensor as following 

no 2w 2(8 -2)w 2 

5°° = — + 4 +0 6, 90 

c z c 4 

9 (H = -^±^ + 0(5), (91) 



(4 7 2 - 2s)w 2 - |] - l q<ij> + 0(6) . (92) 
Then, the Parametrized 2PN Christoffel symbols: 

14 = -^ + 0(5), (93) 

r& = -^ + 2(/3 ' 4 )w ^ +Q(6), (94) 

<r c 4 

n = -^- 2(1+ 4 7Kt + 2(/3 + : )W ^ +0(6), (95) 
c c c 

J% = 2i ^-W)+%^ + 0(5) ) (96) 
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■Ok 



2(1 + 7) 



.^■^+0(5) 



2,(s — 'y jW 6 



Tj 

2~4 W<ij>,fc + Q<ik>,j — Q<jk>,i) + 0(6) . 



(97) 



(98) 



Substituting the 1PN part of Eqs.JHZJ-JSHJ into Eas. l|Mjl and ljfi3|l. we get the parametrized 1PN equations of 
light ray: 



dxip , . 
cn • — r — = —(1 + j)w . 



dt 



cPxip 



(1 + 7) [Vw - 2n(n • Vw)] , 



(99) 
(100) 



which have been shown in Ref.[7] (see Eq.(6.14) and (6.15) of [7]). 

Considering x = x + n(t - to) + x 1P + x 2 p, Eqs.lHU), (|100|l and substituting whole Eas.l|57 )l -l|55 )l into 
Eas. l|64|l and l|65|) we get the parametrized 2PN equations of light ray: 



dx 2 p 2(l + j)w l 



dt 

d Qsqp 

dt 2 z 



(1 + 7). 



1 

2c 



dxu 



-w tn 



dt 

2(1 + 7) 



(/? + 27 + 27 2 - e)w 2 _ eq_ _ mjij) 
c 3 6c 3 2c 3 



k , 4(1 + 7). 



tol' 



2(/3 + 2 7 + 2 7 2 -£) 



-WW , 



2(1 + t) dx 1P „ v < 2(1+t) „ ^1p 
(— - — ■ \7w)n (n- Vw)- 



dt 



dt 



4(g-7 2 -/?+!) 77 k 
2 w(n • Vw)n -q(ij^ k n J n 



3c 2 



i A; 

q,kn n 



V 



e 

6 c 2 



In special case of the solar system, we have g/y\ = 0, therefore the equations become 



(101) 



(102) 



n 



dx 2 p 2(1 + 7)^ 



dt 
d 

~dW z 



(1 + 7) 



1 

2c 



-to t n 



dt 

2(1 + 7) 



(/? + 2 7 + 2 7 2 - e)w 2 



sq_ 

6c 3 



W(j t k)n l n J n + 



* , 4(1+7) 



2(/3 + 27 + 2 7 2 -e) 



-ww , 



2(1 + 7) rfx 1P „ , < 2(1+7) , „ .drip 
( ^ • Vw)n (n- Vmi) — 



4(e- 7 2 -/3 + 1) 



«;(n • Vw)n l 



dt 

^(-v<zK + ^ 



(103) 



(104) 



When /3 = 7 = e = 7j = l, Egs. (E3) , (ESI, fTH^ and lfT0l)l return to Eqs.® , <H2J, <H3J and dHHJ respectively. 
Maybe in our parametrized 2PN equations of light ray, four parameters are not sufficient. But if we compare 
with Ref . , these four parameters are main parameters in the 2PN equations of light ray. 



5 Conclusion Remarks 

1. We have extended DSX scheme to 2PN level approximation, by means of which we can discuss the 2PN 
contribution on the trajectory of light ray. Our extension keep the main advantage in DSX scheme (the linear 
equations of the potential and vector potential, and the linear relationship between w, Wi and W, W a ). In the 
case of the solar system, we obtain a special solution under the PN gauge. But we should point out that if we 
want to discuss the 2PN contribution on the motion of bodies (not light ray), we have to extend the metric 
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further more (goo ~ 0(8) and qoi ~ 0(7)) which someone might finish in future. Considering the precision of 
ASTROD, they will need the theoretical precision up to the 2PN level for the motion of bodies. 

2. We have discussed 2PN equations of light ray based on the extended DSX scheme by means of the 
iterative method. If we have the higher order metric, we also can obtain the higher order equations of light ray 
with a similar iterative method in this paper. Therefore in some meaning the way shown in the section 3 is 
general. The special case for equations of light ray propagating in the solar system(qy = ^q) is also considered. 

3. We also use a parametrized 2PN metric tensor (with four parameters) and its corresponding Christoffel 
symbols to substitute into the null geodetic equations and obtain parametrized 2PN equations of light ray. 
Although we have not induced four parameters strictly, but it agree with the main parameter in |25j . Also if 
2PN part of gij is neglected, the metric tensor return to the one in ^5] which has been deduced strictly. 

Anyway, we hope our work to be useful for the measurement of future space missions. 
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